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Abstract 

A painting consists of objects which are arranged in specific ways. The art of painting is drawing 
the objects, which can be considered as known trends, in an expressive manner. Detrended methods 
are suitable for characterizing the artistic works of the painter by eliminating trends. It means 
that we study the paintings, regardless of its apparent purpose, as a stochastic process. We apply 
multifractal detrended fluctuation analysis to characterize the statistical properties of Mona Lisa, as 
an instance, to exhibit the fractality of the painting. Our results show that Mona Lisa is long range 
correlated and almost behaves similar in various scales. 
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1 Introduction 

Art is a collection of objects created with the intention of transmitting emotions and/or ideas. An 
object can be characterized by intentions of its creator, regardless of its apparent purpose. In this 
sense, art is described as a deliberate process of arrangement by an agent. Art stimulates an individual's 
thoughts, emotions, beliefs, or ideas through the senses. It is also an expression of an idea and it can take 
many different forms and serve many different purposes. Mathematics and art are examples of human 



motivations to understand reality [TJ [2 O |H [6] . The connection between mathematics and art goes 
back to thousands of years. Patterns, symmetries, proportions and transformations are fundamental 
concepts common to both disciplines. One of the connections between mathematics and art is that some 
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known as artists have needed to develop or use mathematical thinking to carry out their artistic insight. 
Thus we are forced to use interpretative techniques in order to search for them. In particular, there 
are several obvious relations between art and mathematics such as music [3j El El [Tj 8 , 9 J and painting 

mm ma mi ei Banana]. 

Painting has three principal parts: drawing, proportion and coloring. Drawing is outlines and con- 
tours contained in a painting. Proportion is these outlines and contours positioned in proportion in their 
places. Coloring is giving colors to things. One can obviously use mathematical geometry to analyze 
paintings, figurative or abstract, in terms of shapes such as points and lines, circles and triangles. One 
may think of the application of the perspective theory to figurative painting, or the use of some concepts 
like fractals for comprehension of abstract paintings. An important example is the analysis of Jackson 
Pollock's drip paintings in terms of fractal geometry [1]. 

Here, we are interested to study the fractal nature of painting from another point of view which in 
general, could be applicable to non-stationary series [3]. Among various paintings we focus on charac- 
terizing the complexity of color signal of "Mona Lisa" through computation of signal parameters and 
scaling exponents, which quantifies correlation exponents and multifractality of the signal. 

Mona Lisa, or La Gioconda (La Joconde), is a 16th century oil painting on poplar wood by Leonardo 
da Vinci, and is the most famous painting in the world. This painting is a half-length portrait which 
depicts a woman whose gaze meets the viewer's with an enigmatic expression (Fig. [1]). Because of 
non-stationary nature of color signal series, and due to finiteness of available data sample, we should 
apply methods which are insensitive to non-stationarities, like trends. In order to separate trends from 
correlations we need to eliminate trends in our color data. Several methods are used effectively for 
this purpose: detrended fluctuation analysis (DFA) [15], rescaled range analysis (R/S) [16] and wavelet 
techniques (WT) [17j . However, using increment series is customary to make a stationary series from a 
non-stationary ones |18j . 

We use MF-DFA method for analysis and eliminating trends from data set. This method is the 
modified version of DFA method to detect multifractal properties of time series. DFA method introduced 
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Figure 1: Mona Lisa: oil painting on poplar wood by Leonardo da Vinci. 

by Peng et al. [15] has become a widely used technique for the determination of monofractal scaling 
properties and the detection of long-range correlations in noisy, non-stationary time series [19, 20, 2Hl22j. 
It has successfully been applied to diverse fields [TilllllllllSllMllSTllllMlllSllSniE]- One reason to 
employ DFA method is to avoid spurious detection of correlations that are artefacts of non-stationarity 
series. The focus of the present paper is on the fractality nature of color series obtained from Mona Lisa. 
To construct the series, we calculate the standard color values of each pixel of the picture successively 
for each row from up to down, continuously. In particular, Fig. [T] shows Mona Lisa and Fig. [5] shows its 
gray, red, green and blue color styles for a 202 x 300 pixels sample. The color fluctuation graphs of gray 
and red styles, are also depicted in Fig. [31 

The paper is organized as follows: In Sec. [2[ we describe MF-DFA methods in detail and show 
that scaling exponents determined by MF-DFA method are identical to those obtained by standard 
multifractal formalism based on partition functions. In Sec. [3j in analysis of color series of Mona Lisa 
we also examine the multifractality in color data. Section [4] closes with a conclusion. 
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2 Multifractal Detrended Fluctuation Analysis 

The simplest type of multifractal analysis is based upon standard partition function multifractal formal- 
ism, which has been developed for multifractal characterization of normalized, stationary measurements 
[3H [331 1311 135] ■ Unfortunately, this standard formalism does not give us correct results for non-stationary 
time series that are affected by trends or those which cannot be normalized. MF-DFA is based on iden- 
tification of scaling of gth-order moment depending on signal length, and this is a generalization of 
standard DFA method in which q = 2. Moreover one should find correct scaling behavior of fluctuations, 
from experimental data which are often affected by non-stationary sources, like trends. These have to 
be well distinguished from intrinsic fluctuations of the system. In addition, often in collected data we do 
not know the reasons, or even worse the scales, for underlying trends, and also available record data is 
usually small. So, for a reliable detection of correlations, it is essential to distinguish trends for intrinsic 
fluctuations from collected data. Hurst rescaled-range analysis |16j and other non-detrending methods 
work well when records are long and do not involve trends, otherwise they might give wrong results. 
DFA is a well established method for determining scaling behavior of noisy data where the data include 
trends and their origin and shape are unknown [TS1 HSJ [351 133 H3j • 

Modified multifractal DFA procedure consists of five steps. The first three steps are essentially 
identical to conventional DFA procedure (see e.g. [HI [HI [2Ql EH [22] ) ■ Suppose that Xk is a series of 
length N , and it is of compact support, i.e. Xk — for an insignificant fraction of the values only. 
• Step 1: Determine the profile 



Subtraction of the mean (x) from Xk is not compulsory, since it would be eliminated by later detrending 
in third step. 

• Step 2: Divide profile Y(i) into N s = int(N/s) non-overlapping segments of equal lengths s. Since 
length N of series is often not a multiple of considered time scale s, a short part at the end of profile may 
remain. In order not to disregard this part of the series, same procedure should be repeated starting 
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Figure 3: Typical color value of color series for gray (Up) and red (Down) style. 
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from the opposite end. 

• Step 3: Calculate local trend for each of N s segments by a least-square fit of the series. Then determine 
the variance 

F 2 ( S , V ) = \± {Y[{v l)s + i] - y„(i)} 2 , (2) 
i=i 

for each segment v = 1, . . . ,N S . Where, y u {i) is fitted polynomial in segment v. Linear, quadratic, cubic, 
or higher order polynomials can be used in fitting procedure (conventionally called DFA1, DFA2, . . ., 
DFAto) [151 126] . In (MF-)DFAm, trend of order m in profile (and equivalently, order m — 1 in original 
series) are eliminated. Thus a comparison of results for different orders of DFA allows one to estimate 
the type of the polynomial trend in the series [20J, [21] . 

• Step 4'- Average over all segments to obtain q-th order fluctuation function, defined by: 

where, in general, variable q can take any real value except zero. For q = 2, standard DFA procedure is 
retrieved. Generally we are interested to know how generalized q dependent fluctuation functions F q (s) 
depend on time scale s for different values of q. Hence, we must repeat steps 2, 3 and 4 for several scales 
s. It is apparent that F q (s) will increase with increasing s. Of course, F q (s) depends on DFA order m. 
By construction, F q (s) is only defined for s > m + 2. 

• Step 5: Determine scaling behavior of fluctuation functions by analyzing log- log plots of F q (s) versus 
s for each value of q. If series Xj are long-range power law correlated, then F q (s), for large values of s, 
increases as a power- law i.e., 

F q (s)~s h M. (4) 

In general, exponent h(q) may depend on q. For stationary series such as fractional Gaussian noise 
(fGn), Y(i) in Eq. [T]will have a fractional Brownian motion (fBm) signal, so, < h(q = 2) < 1.0. The 
exponent h(2) is identical to well known Hurst exponent H [15l [T9l [32] . Also, for non-stationary signals, 
such as fBm noise, Y(i) in Eq. [T]will be a sum of fBm signal, so corresponding scaling exponent of F q (s) 
is identified by h(q — 2) > 1.0 [T51 138j . For monofractal series, h(q) is independent of q, since scaling 
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behavior of variance F 2 (s, v) is identical for all segments v, and averaging procedure in Eq. ([3]) will just 
give a same scaling behavior for all values of q. If we consider positive values of q, the segments v with 
large variance F 2 (s,v) (i.e. large deviation from the corresponding fit) will dominate average F q (s). 
Thus, for positive values of q, h(q) describes scaling behavior of segments with large fluctuations. For 
negative values of q, segments v with small variance F 2 (s,v) will dominate average F q (s). Hence, for 
negative values of q, h(q) describes scaling behavior of segments with small fluctuations. 

For a stationary and normalized scries, multifractal scaling exponent h(q) defined in Eq. (JJ) is directly 
related to scaling exponent r(q) defined by standard partition function based on multifractal formalism 
with the following analytical relation (i.e. see [32]) 

T(q)=qh(q)-1. (5) 

Thus, we see that h(q) defined in Eq. ([4]) for MF-DFA is directly related to classical multifractal scaling 
exponent r(q) and generalized multifractal dimension [32 

D(q) = ^ = (6) 

that is used instead of r(q) in some papers. In this case, while h(q) is independent of q for a monofractal 
series, D(q) depends on q. Another way to characterize a multifractal series is looking to singularity 
spectrum /(a), which is related to r(q) via a Legendre transform [3"2l l3~i] 

a = t' (q) and f(a)—qa — T(q). (7) 

Here, a is singularity strength or Holder exponent, while f(a) denotes dimension of a subset of series 
that is characterized by a. Using Eq. ([5]), we can directly relate a and f(a) to h(q) 

a = h(q) + qh' (q) and f(a)=q[a-h(q)] + l. (8) 

Holder exponent denotes monofractality, while in multifractal case, different parts of structure are char- 
acterized by different values of a, leading to existence of spectrum f(a). 



8 



10' 



10 s 

S 



10' 



10 z 10 J 10' 

s 



Figure 4: The log-log plot F(s) versus s for q — 2.0 after double profiling for gray style, before F-DFA 
(Left) and After F-DFA (Right). 

2.1 Fourier-Detrended Fluctuation Analysis 

In some cases, there exists one or more crossover length scales si in graph F(s), separating regimes 
with different scaling exponents [20 . In these cases, obtaining the scaling behavior is more complicated 
and different scaling exponents are required for different parts of the series 22 . Therefore, one needs 
a multitude of scaling exponents in various scales for a full description of scaling behavior. A crossover 
usually can arise from a change in correlation properties of the signal at different scales, or can often arise 
from trends in data. In addition, inappropriate detrcnding can create additional trends to signal which 
appear as crossovers in F(s). The extra crossovers are not actual length scales of the original signal. 
In fact, these crossovers are due to the unsuitable detrending process. Fourier-Detrend Fluctuation 
Analysis (F-DFA) can be applied to remove crossovers such as sinusoidal trends. The F-DFA is a 
modified approach for the analysis of low frequency trends [39l [40], [41] . 

We transform data record to Fourier space in order to remove trends having a low frequency periodic 
behavior. Then we truncate the first few coefficients of the Fourier expansion and inverse Fourier 
transform of the series. After removing the sinusoidal trends we can obtain the fluctuation exponent 
using the direct calculation of the MF-DFA. If truncation numbers are sufficient, The crossovers due to 
sinusoidal trends in the log- log plot of F q (s) versus s disappear. 
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3 Analysis of Color series 

Most paintings are non-stationary due to the fact that all information is in whole of their structures. 
It means if we cut a part of a painting, we can not complete the remainder of it. Although there 
are some fractal paintings which all information is in each part of them. As mentioned in Sec. [H a 
spurious of correlations may be detected using detrending methods in non-stationarity series. While, 
direct calculation of correlation behavior, spectral density exponent, fractal dimensions, etc., do not give 
us reliable results. 

It seems that most of non-stationary properties in a painting are known trends such as objects. In 
fact, a painting consists of some objects or components have been set in various parts of it. As an 
example, a portrait has different parts such as eyes, nose, mouth, cheeks, eyebrows, etc. which can 
be considered as known trends. Also, light and shadow can be considered as trends. Since the art of 
painting is drawing the objects in an expressive manner, we intend to detrend data to characterize the 
artistic work of a painter regardless of the existence of specific objects (trends) in his or her paintings. 

To construct color series, we need to relate a standard numerical quantity to each pixel. The basic 
colors are Red, Green and Blue (RGB), which every color is a specific combination of these basic colors. 
The standard value of each basic color is in the range of — 255. Therefore, we consider each color in 
our analysis and also we cover the gray style, which is a certain combination of these three basic colors. 
In the first step, it can be checked out that, all color series are non-stationary. One can verify non- 
stationarity property experimentally by measuring stability of their average and variance in a moving 
window for example using scale s. In fact, there are some specific lengths in painting that have important 
effects in statistical parameters. The picture's width and height and color intensities are the first scaling 
parameters in painting. Moreover, the body and face detail structures are important scales that painter 
usually regards them. Let us determine whether the data set has a sinusoidal trend or not. According to 
MF-DFA1 method, Generalized Hurst exponents h(q) in Eq. (fj| can be found by analyzing log-log plots 
of F q (s) versus s for each q (Fig. 2}. Our investigation shows that there is one crossover length scale si 
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in log- log plots of Fq(s) versus s for every q's. To cancel the sinusoidal trends in MF-DFAl, we apply 
F-DFA method for color data. For eliminating the crossover scales, we need to remove only one term 
of the Fourier expansion. Then, by inverse Fourier Transformation, the noise without sinusoidal trend 
is extracted. Hurst exponent is between < H < 1. However, MF-DFA method can only determine 
positive generalized Hurst exponents, in order to refine the analysis near the fGn/fBm boundary or 
strongly anti-correlated signals when it is close to zero. 

The simplest way to analyze such data is to integrate series before MF-DFA procedure. Hence, 
we replace the single summation in Eq. [T] which is describing the profile from the original data, by a 
double summation using signal summation conversion method (SSC) 38, 42l|43l[3j. After using SSC 
method, fGn switch to fBm and fBm switch to sum-fBm. In this case the relation between the new 
exponent, h(q = 2), and H is H = h(q = 2) — 1 [38l|42l|43] (recently Movahed et al. have proven the 
relation between derived exponent from double profile of series in DFA method and h(q — 2) exponent 
in Ref. [H]). We find h(q = 2) = 1.74 ± 0.01 for gray style color series using SSC method. Therefore, 
Hurst exponent equals to H = h(q = 2) — 1 = 0.74. 

The results of MF-DFAl method for color signal are shown in Fig. 03 which show that these color series 
can be considered approximately as a monofractal process which is indicated by weak q dependence of the 
exponents h(q — 2) and r(q) [43] . The q dependence of multifractal scaling exponent r(q) approximately 
has a linear dependence to q with equal slopes as 1.63, 1.59, 1.70, and 1.66 for gray, red, green, and blue 
styles, respectively. Table [1] shows the obtained quantities using MF-DFAl method. Figure [5] shows the 
width of singularity spectrum (/(a)) for gray style with Aa ~ 0.15 (Aa = a(q max ) — a(q m i n )). Its value 
indicates that the power of multifractality of the color series is weak |44) . 

4 Conclusions 

Using detrended methods for analyzing paintings, we can study the artistic manner of the painter, 
regardless of subject of the artworks. This is due to the fact that a painting consists of objects which 
can be considered as trends. The non-stationary property of paintings can be related to the objects, 
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Table 1: Values of h(q = 2), r(q = 2) exponents and width of singularity spectrum, /(a), Aa for q = 2.0 
of various styles obtained by MF-DFA1. 



Style 


h(2) 


r(2) 


Aa 


Gray 


1.74 ±0.01 


2.48 ±0.01 


0.15 


Red 


1.74 ±0.01 


2.48 ±0.01 


0.11 


Green 


1.72 ±0.01 


2.44 ±0.01 


0.12 


Blue 


1.74 ±0.01 


2.48 ±0.01 


0.12 



light and shadow, etc. Indeed by detrending methods, we eliminate trends from paintings and make the 
signal stationary. It means, we reduce the effects of objects, light and shadow, etc. from paintings. We 
have used multifractal detrended fluctuation analysis to detrend data to characterize the most famous 
artistic work of Leonardo da Vinci and shown the fractality of this painting. Our results using MF-DFA 
show that color series of Mona Lisa almost has the same behavior in various scales for all studied color 
styles (red, green, blue and gray). 
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